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^ ■ We present a detailed chemical abundance study of seven giant stars in M68 

^ ■ including six red giants and one post-AGB star. We find significant differences 

\Q • in the gravities determined using photometry and those obtained from ionization 

^ ■ balance, which suggests that non-LTE effects are important for these low-gravity, 

Q\ • metal-poor stars. We adopt an iron abundance using photometric gravities and 

^ ' Fe II lines to minimize those effects, finding [Fe/H] = —2.16 ± 0.02 (cr = 0.04). 

O ■ For element-to-iron ratios, we rely on neutral lines vs. Fe I and ionized lines vs. 

1^ • Fe II (except for [0/Fe]) to also minimize non-LTE effects. 

We find variations in the abundances of sodium among the program stars. 

2 ' However, there is no correlation (or anti-correlation) with the oxygen abundances. 

^ ■ Further, the post-AGB star has a normal (low) abundance of sodium. Both 
■ 

^ ' of these facts add further support to the idea that the variations seen among 

v~j ■ some light elements within individual globular clusters arises from primordial 

■ variations, and not from deep mixing. 



M68, like M15, shows elevated abundances of silicon compared to other glob- 
ular clusters and comparable metallicity field stars. But M68 deviates even more 
in showing a relative underabundance of titanium. We speculate that in M68, 
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titanium is behaving like an iron-peak element rather than its more commonly 
observed adherence to enhancements seen in the "a" elements such as magne- 
sium, silicon, and calcium. We interpret this result as implying that the chemical 
enrichment seen in M68 may have arisen from contributions from supernovae with 
somewhat more massive progenitors than contribute normally to abundances seen 
in other globular clusters. 

The neutron capture elements barium and europium vary among the stars 
in M15 (Sneden et al. 1997), but the [Ba/Eu] is relatively constant, suggesting 
that both elements arise in the same nucleosynthesis events. M68 shares the 
same [Ba/Eu] ratio as the stars in M15, but the average abundance ratio of these 
elements, and lanthanum, arc lower in M68 relative to iron than in M15, implying 
a slightly weaker contribution of r-process nucleosynthesis in M68. 

Subject headings: Galaxy: halo — globular clusters: individual (M68:NGC 4590) 
— stars: abundances 



1. INTRODUCTION 

Globular clusters in our Galaxy do not define a single homogeneous population with a 
single history. In his pioneering work, Zinn (1985) subdivided the globular clusters into the 
"halo" and the "thick disk" groups at [Fe/H] = —0.8. The halo clusters have an essentially 
spherical distribution about the Galactic center and they constitute a pressure supported 
system (a small rotational velocity and a larger velocity dispersion), while the thick disk 
clusters have a highly flattened spatial distribution and constitute a rotational supported 
system (a larger rotational velocity and a smaller velocity dispersion). Searle & Zinn (1978) 
and Lee, Demarque, Sz Zinn (1994) suggested that the inner halo globular clusters exhibit 
a tight horizontal branch (HB) morphology versus [Fe/H] relation, while the outer halo 
globular clusters show the second parameter phenomenon (i.e., a larger scatter in HB type^ 
at a given [Fe/H]). Subsequently, Zinn (1993) subdivided the halo clusters into two groups. 
The "old halo" group obeys the same HB type versus [Fe/H] relationship as the inner halo 
clusters while the "younger halo" group deviates from this relationship by a significant 
amount (sec Figure 1). Zinn (1993) and Da Costa & Armandroff (1995) argued that the 
old and the younger halo groups have different kinematic properties that the old halo group 
has a prograde mean rotation velocity with a smaller velocity dispersion while the younger 



^HB type is defined to be (B— R)/(B-|-V-|-R), where B, V, R are the numbers of blue HB stars, RR Lyrae 
variables, red HB stars, respectively (Lee, Demarque, & Zinn 1990). 
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halo group has a retrograde mean rotation velocity about the Galactic center with a larger 
velocity dispersion. They suggested that the old halo group formed during the collapse that 
led ultimately to the formation of the Galactic disk and the younger halo group were accreted 
later in time. If so, one would expect to see signatures of different chemical enrichment 
history carved in spectra of stars in globular clusters. For example, one might be able to 
investigate the early form of the initial mass function (IMF) by studying abundances of 
stellar mass-sensitive elements (see discussion and references in McWiUiam 1997). 

An alternative perspective on "younger halo" and "old halo" clusters is that of dissolved 
and accreted dwarf galaxies. For example. Freeman (1993) suggested that the massive and 
chemically unusual globular cluster uj Gen might be the remnant nucleus of an accreted dwarf 
galaxy. Lynden-Bell & Lynden-Bell (1995) noted the possible ahgnments of orbital poles of 
some globular clusters such that they might comprise a "spoor" (more commonly referred to 
as a stream) of clusters. While the specific details have not found support, the idea has been 
demonstrated very nicely by Dinescu et al. (2000), who found a clear dynamical relationship 
between the space motions of the young globular cluster Palomar 12 and the Sagittarius 
dwarf galaxy. Indeed, other globular clusters may also be associated with this particular 
accretion event (NGG 5634: Bellazzini, Ferraro, & Ibata 2002; NGG 4147: Bellazzini et al. 
2003; Palomar 2: Majewski et al. 2004). Yoon & Lee (2002) have similarly speculated on a 
dynamical relationship between some metal-poor globular clusters, including M68. Yoon & 
Lee (2002) suggested that M68 was a member of a sateUite galaxy and accreted later in time 
to our Galaxy, resulting in a planar motion of several metal-poor clusters including M15 and 
M92. 

Dynamical evolution in the Galaxy can eventually dissolve streams, making their de- 
tection difficult. However, an interesting alternative exists, called "chemical tagging". As 
Freeman & Bland-Hawthorn (2002) discussed, stars born in galaxies whose star formation 
histories differ from those that have created the bulk of the Galaxy's stars may still be dis- 
cernible in unusual clcment-to-iron ratios. Indeed, Cohen (2004) has found a compelling link 
between Palomar 12 and the Sagittarius dwarf. Detailed chemical abundances of globular 
clusters may yet become the principle means of identifying historical links between stars and 
clusters that are now widely dispersed in our Galaxy. 

M68 (NGG 4590) is part of the younger halo, according to its HB type with respect to its 
metallicity, and is located 10 kpc from the Galactic center. Dinescu, Girard, & van Altena 
(1999) have found that the cluster's Galactic orbit carries it rather far from the Galactic 
center, perhaps as far as 30 kpc. As shown in Figure 1, globular clusters with [Fe/H] < —2.0 
have HB types that are greater than 0.6 with the exception of M68. If age is the second 
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parameter^, M68 represents the low metallicity tail of the younger halo globular clusters. The 
absolute age of M68 does appears to be slightly younger than those of the oldest globular 
clusters in our Galaxy. Rosenberg et al. (1999) argued that M68 is coeval to or slightly 
younger than those of the oldest halo globular clusters while VandenBerg (2000) discussed 
that the age of M68 is less than that of M92, one of the oldest globular clusters in our Galaxy, 
by ^ 15%. 

The previous metallicity estimates of M68 as follows. Zinn & West (1984) and Zinn 
(1985) adopted [Fe/H] = -2.09. Gratton & Ortolani (1989) observed two stars in M68 
with one being in common with this study. They derived elemental abundances for 11 
species including iron, finding [Fe/H] = —1.92. However, due to the low resolving power 
{R — 15,000) of their spectra, their equivalent width measurements were vulnerable to line 
blending. Minniti et al. (1993) observed two stars in M68, one being in common with 
Gratton & Ortolani (1989) and this study, finding [Fe/H] = -2.17. Minniti et al. (1996) also 
discussed oxygen and sodium abundances of the cluster. Finally, Rutledge (1997) observed 
the infrared Ga II triplet for 19 stars and obtained [Fe/H] = —2.11 ± 0.03 using the Zinn & 
West (1984) [Fe/H] scale. 

In this paper, we explore the detailed elemental abundances for seven giant stars in M68. 

One of these stars is a probable asymptotic giant branch (AGB) stars and the remaining six 
are red giant branch (RGB) stars. This study is tied directly to that of Lee & Garney (2002) 
with the same instrument setups and analysis methods. 

2. OBSERVATIONS AND DATA REDUCTION 

The observations were carried out from 4 to 7 May 1996. We selected our program 
RGB stars from the BV photometry of Walker (1994). The positions of our target giant 
stars on the color-magnitude diagram along with bright stars in M68 are shown in Figure 2. 
In Table 1, we provide identifications (Alcaino 1977; Harris 1975; Walker 1994), position 
(Gutri et al. 2000), V magnitudes, (B — V) colors (Walker 1994) and K magnitudes (Frogel, 
Persson, & Gohen 1983; Gutri et al. 2000) of our target stars. Please note that the 2MASS K 
magnitudes have been converted to the GIT system (Gutri et al. 2000) are in good agreement 
with those of Frogel et al. (1983). We obtained high signal- to- noise ratio (S/N > 90 per pixel) 
echelle spectra using the GTIO 4-meter telescope and its Gassegrain echelle spectrograph. 
The Tek 2048 x 2048 GGD, 31.6 lines/mm echelle grating, long red camera, and G181 



^Thc "second parameter" is the disputed variable that, in addition to mean metallicity, determines the 
HB type of a cluster. 
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cross-disperser were employed for our observations. The slit width was 150 //m, or about 
1.0 arcsec, that projected to 2.0 pixels and which yielded an effective resolving power R — 
28,000. Each spectrum had complete spectral coverage from 5420 to 7840 A. All program 
star observations were accompanied by flat lamp, Th-Ar lamp, and bias frames. 

The raw data frames were trimmed, bias-corrected, and flat-fielded using the IRAF^ 
ARED and CCDRED packages. The scattered light was also subtracted using the AP- 
SCATTER task in ECHELLE package. The echelle apertures were then extracted to form 
1-d spectra, which were continuum-fitted and normalized, and a wavelength solution was 
applied following the standard IRAF echelle reduction routines. 

Equivalent widths were measured mainly by the direct integration of each line profile 
using the SPLOT task in IRAF ECHELLE package. We estimate our measurement error in 
equivalent width to be ±2 mA from the size of noise features in the spectra and our ability 
to determine the proper continuum level. The equivalent widths for our program stars are 
hsted in Table 2. 

Gratton & Ortolani (1989) and Minniti et al. (1993) obtained a spectrum of one star in 
common with our program stars. They used the identifications from Harris (1975) while we 
have used those of Walker (1994). Comparing the two shows that the star 1-260 from Harris 
(1975) is the same as the star 160 from Walker (1994). Figure 3 compares equivalent width 
measurements of our work with those measured by Gratton & Ortolani (1989; crosses) and 
Minniti et al. (1993; open circles). The agreement with the latter study is quite good, with 
a mean difference of 2.3 ± 1.5 mA (in the sense of their study minus ours). The instrumental 
resolving powers for the two studies were very similar (27.000 vs. 28,000). However, the 
lower resolving power of the observations reported by Gratton & Ortolani (1989), about 
15,000, appears to have led to systematically larger equivalent widths, as Figure 3 reveals. 
The mean difference is 13.8 ± 2.7 mA. 

3. ANALYSIS 

In our elemental abundance analysis, we use the usual spectroscopic notations that 
[A/B] = \og{NA/NB)star ' log{NA/NB)Q, and that \ogn{A) = logiNA/Nn) + 12.00 for each 
element. For the absolute solar iron abundance, we adopt logn(Fe) = 7.52 following the 



•^IRAF (Image Reduction and Analysis Facility) is distributed by the National Optical Astronomy Obser- 
vatory, which is operated by the Association of Universities for Research in Astronomy, Inc., under contract 
with the National Science Foundation. 
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discussion of Sneden et al. (1991). 

3.1. Line Selection and Oscillator Strengths 

For our line selection, laboratory oscillator strengths were adopted whenever possible, 
with supplemental solar oscillator strength values. In addition to oscillator strengths, taking 
into account the damping broadening due to the van der Waals force, we adopted the Unsold 
approximation with no enhancement. 

The abundance analysis depends mainly on the reliability of the oscillator strength 
values of the Fe I and Fe II lines, since not only the metaUicity scale but also the stellar 
parameters, including the spectroscopic temperature, surface gravity, and microturbulent 
velocity, will be determined using these hues. As discussed by Lee & Carney (2002), we 
mainly relied upon the extensive laboratory oscillator strength measurements by the Oxford 
group (Blackwell et al. 1982b, 1982c, 1986a). We also used oscillator strength values mea- 
sured by O'Brian et al. (1991) and the Hannover group (Bard, Kock, & Kock 1991; Bard & 
Kock 1994). In our iron abundance analysis, wc consider the Oxford group's measurements 
(the absorption method) as the "primary" oscillator strengths and oscillator strength mea- 
surements that relied on emission methods (O'Brian et al. 1991; Bard, Kock, & Kock 1991; 
Bard & Kock 1994) as "secondary". Therefore, the oscillator strengths by O'Brian et al. and 
the Hannover group were scaled with respect to those by the Oxford group as a function of 
excitation potential (de Almeida 2000, private communication), 

hggf = log (7/(05) -0.017, 

\oggf = log(7/(if91)- 0.015 - 0.009X, 

loggf = log5/(i/94)- 0.027 - 0.009X, (1) 

where the excitation potential x is given in electron volts. Blackwell, Smith, & Lynas-Gray 
(1995) also pointed out that there appears to exist a shght gradient in the excitation potential 
between oscillator strengths by the Oxford group and those by the Hannover group, with 
loggf (Oxford) — loggf (Hannover) — 0.021 — 0.006%. 

For neutral titanium lines, we relied on the laboratory measurements by the Oxford 
group (Blackwell et al. 1982a, 1983, 1986b). It should be noted that the original Oxford gf- 
values have been increased by +0.056 dex following Grevesse, Blackwell, & Petford (1989). 
They discussed that the Oxford 5f/-values relied on the inaccurate lifetime measurements 
and the absolute /-values should be revised based on the new measurements. 

Hyperfine splitting (HFS) components must be considered in the barium abundance 
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analysis because Ba II lines are usually very strong even in metal-poor stars and the desat- 
uration effects due to HFS components become evident (see for example, McWilliam 1998). 
We adopted the Ba II HFS components and oscillator strengths of Sneden et al. (1997). 
We also perform HFS treatment for scandium, manganese (Prochaska & McWilliam 2000) 
and copper (Kurucz 1993). For the copper HFS analysis, we adopt the solar Cu isotopic 
ratio, 69% ^^Cu and 31% ^^Cu, following the discussions given by Smith et al. (2000) and 
Simmerer et al. (2003). The equivalent widths of Mn I A 6021.79 A and Cu I A 5782.13 
A lines in our program stars are weak and the Mn and Cu abundance differences between 
HFS treatment and non-HFS treatment are no larger than 0.02 dex. We list our source of 
oscillator strengths for each element in Table 3. 



3.2. Stellar Pctrameters and Model Atmospheres 

Having good stellar parameters, such as the effective temperature and the surface grav- 
ity, is critical for any stellar abundance study, since the absolute or the relative elemental 
abundance scale will depend on the input stellar parameters. For our analysis, we rely 
on spectroscopic temperatures and photometric surface gravities, following the method de- 
scribed in Lee & Carney (2002; see also Ivans et al. 2001, Kraft & Ivans 2003, and Sneden et 
al. 2004). It has been suspected by others that the traditional spectroscopic surface gravity 
determination method which requires the same elemental abundances derived from neutral 
and singly ionized lines (preferentially Fe I and Fe II lines) suffers from non-local thermody- 
namic equilibrium (NLTE) effects (see, for example, Nissen et al. 1997; AUende Prieto et al. 
1999). Since metal-poor stars have much weaker metal-absorption in the ultraviolet, more 
non-local UV flux can penetrate from the deeper layers. This flux is vital in determining the 
ionization equilibrium of the atoms, resulting in deviations from local thermodynamic equi- 
librium (LTE). Nissen et al. (1997) claimed that surface gravities of metal-poor dwarfs and 
subgiants derived from the spectroscopic method, which demands that Fe I and Fe II lines 
should provide the same iron abundance, arc a factor of two or three (Aloggf ~ 0.3 - 0.5) 
smaller than those from the Hipparcos parallaxes. Allcndc Pricto et al. (1999) also claimed 
that spectroscopic gravities and those from the Hipparcos parallaxes are in good agreement 
for stars in the metallicity range —1.0 < [Fe/H] < -1-0.3, while large discrepancies can be 
found for stars with metallicities below [Fe/H] = —1.0, in the sense that the spectroscopic 
method provides lower surface gravities. Therefore, we rely on photometric gravities for our 
abundance analysis. 

The initial estimates of the temperature of program stars were estimated using BVK 
photometry of our program stars (Cutri et al. 2000; Frogel, Persson, & Cohen 1983; Walker 
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1994) and the empirical color-temperature relations given by Alonso, Arribas, & Martinez- 
Roger (1999). Since their relation depends slightly on the metallicity, we adopted [Fe/H] 
= —2.1 for M68 (Harris 1996). To estimate the dereddened color, we adopt E{B — V) = 
0.07 (Walker 1994). Note that Alonso et al. (1999) employed the Carlos Sanchez Telescope 
(TCS) system for their infrared color-temperature relations and we used the relation given 
by Alonso, Arribas, & Martinez- Roger (1998) to convert the CIT system to the TCS system 
for K magnitudes hsted in Table 1. To derive photometric surface gravity in relation to that 
of the Sun, we use loggQ — A^AA in cgs units, Mboi,© = 4.74 mag, and T^s^q — 5777 K for 
the Sun (Livingston 1999) and we assume the stellar masses for all the red giants stars in 
this analysis to be M = 0.8 Mq. We use the empirical relation given by Alonso et al. (1999) 
to estimate the bolometric correction and we adopt (m — M)o = 14.97 mag for the cluster 
(Harris 1996). 

With initial photometric temperature and surface gravity estimates, 72-depth plane- 
parallel LTE model atmospheres were computed using the program ATLAS9, written and 
supplied by Dr. R. L. Kurucz. Assuming that all of the cluster's stars would prove to be 
metal-poor, the model atmospheres were computed using opacity distribution functions and 
abundances with enhanced abundances of all the "a" elements (O, Ne, Mg, Si, S, Ar, Ca, 
and Ti) by 0.4 dex. The "a" element enhancements are important since several of these 
elements are quite abundant and are major electron donors to the H~ opacity. During our 
model computation, convective overshoot was turned off. 

The abundance analysis was performed using the current version (2002) of the LTE 
line analysis program MOOG (Sneden 1973). Adopting the photometric temperature and 
surface gravity as our initial values, we began by restricting the analysis to those Fe I lines 
with log(WA/A) < —5.2 (i.e., for the linear part of the curve of growth), and comparing the 
abundances as a function of excitation potential. New model atmospheres were computed 
with a slightly different effective temperature until the slope of the log n(Fe I) versus excita- 
tion potential relation was zero to within the uncertainties. The stronger Fe I lines were then 
added and the microturbulent velocity t'turb altered until the log n(Fe I) versus log(WA/A) 
relation had no discernible slope. 

Table 4 shows our temperature and surface gravity of program stars. In the fifth column 
of the Table, we also show the temperature of five stars in common given by Frogel et al. 
(1983). Please note that Gratton & Ortolani relied on temperatures given by Frogel et al. 
(1983), which are 159 ± 27 K higher than our spectroscopic temperatures. The temperature 
difference between those derived from {B — V) colors (column 2) and spectroscopic temper- 
atures of this study is 72 ± 15 K (7 stars), in the sense that our spectroscopic temperature 
is low. The temperature difference between those derived from (V — K) colors using K 
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magnitudes of Progel et al. (column 3) and spectroscopic temperatures of this study is 54 ± 
17 K (5 stars) and those derived from {V — K) colors using 2MASS K magnitudes (column 
4) and spectroscopic temperatures of this study is 69 ± 14 K (7 stars). The discrepancy 
between photometric and spectroscopic temperatures from [V — K) colors is slightly smaller 
than that from [B — V) colors. 

4. RESULTS 

4.1. Elemental Abundances and Error Analysis 

In Tables 5, we present the elemental abundances of our program stars using photometric 
surface gravities and spectroscopic temperatures. The [el/Fe] ratios for neutral elements are 
estimated from [el/H] and [Fe I/H] ratios, with the exception of oxygen. The [el/Fe] for 
singly ionized elements (Ti II, Ba II, La II, and Eu II) and oxygen are estimated from [el/H] 
and [Fe II/H] ratios (see, for example, Ivans et al. 2001, Kraft &: Ivans 2003, Sneden et al. 
2004). The internal uncertainty quoted is for a single line and, therefore, that of each clement 
is given by a/^/n, where a is the uncertainty per line and n is the number of absorption 
lines used for each element. Systematic errors, such as in adopted gf values as a function 
of excitation potential, which could lead to systematically erroneous temperature estimates, 
are not included. The last two columns of the Table show the mean values of each element 
of the cluster with and without the star 117 (ZNG2), an ultraviolet-bright, post-AGB star 
according to Zinn, Newell, & Gibson (1972). We adopt iron abundances based only on the 
Fe II lines for our program stars since the Fe II abundance is thought to be less sensitive 
to NLTE conditions (see for example, Thevenin & Idiart 1999, Kraft & Ivans 2003). In 
Figure 4, we show abundances of all elements measured in this study against Teg, showing 
no discernible gradient in elemental abundances with Tcff. The mean [Fe/H] of —2.16 dex 
for our seven stars is measured with a small internal uncertainty of ±0.02 (cr = 0.04). 

For comparison, we also show the elemental abundances of our program stars using 
the traditional spectroscopic surface gravities (see column 7 of Table 4). In the Table, the 
[el/Fe] ratios are derived from [el/H] and mean [Fe/H] ratios. In Table 7, we show differences 
in elemental abundances using photometric and spectroscopic gravities. The [Fe/H] ratio 
using spectroscopic gravities is 0.25 dex lower than that using photometric gravities. For 
other elements, however, the elemental abundances from two different methods are in good 
agreement to within 0.10 dex. Therefore, these small differences suggest that we are still 
able, in principle, to compare our elemental abundances with other results using the different 
surface gravity determination method. Nonetheless, we reiterate that we choose to employ 
photometric gravities and will discuss only the results from Table 5. 
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In Table 8, we show estimated errors resulting from uncertainties in the input model 
atmosphere ST^s = ± 80 K, Slogg = ±0.3 and 5fturb — ±0.2 km s~^, which are appropriate 
for our analysis. The Fe, Si, Ba, and La abundances are sensitive to Teg resulting in |(5[el/Fe 
(or H for Fe)]/5Tcfj(80 K)| ^ 0.08 - 0.12 dex. The iron abundance from Fe II lines is sensitive 
to surface gravity, \6[Fe II/H]/51og5'(0.3 dex)| ^ 0.10 dex. Since our program stars are metal- 
poor, absorption lines are usually weak as listed in Table 2. Therefore, our derived elemental 
abundances are less sensitive to the microturbulent velocity. The barium abundance is the 
most sensitive to the microturbulent velocity resulting in |5[Ba/Fe]/5fturb(0.2 km s~^)| 
0.09 dex. 



4.2. Compcirisons with previous results 

As mentioned above, Gratton & Ortolani (1989) and Minniti et al. (1993, 1996) derived 
abundances for star 160 (1-260 from Harris 1975). In Table 9, we compare our stellar param- 
eters and elemental abundances of this star with those of Gratton & Ortolani (1989) and 
Minniti et al. (1993, 1996). 

Gratton & Ortolani (1989) relied on the results from Frogel et al. (1983) for the estimated 
temperature and surface gravity. Table 4 shows that this results in a Tes estimate 230 K 
hotter than our spectroscopic measurement, and a slightly lower microturbulent velocity (by 
0.2 km s~^) as well. Our results suggest that a difference in microturbulent velocity of 0.2 
km s~^ does not alter the derived iron abundances and elemcnt-to-iron ratios by more than 
about 0.1 dex. As discussed above (see also Figure 3), the equivalent widths measured by 
Gratton & Ortolani are about 14 mA larger than our measurements. This was also noticed 
by Minniti et al. (1993), who found a difference in equivalent width of 10.2 ± 9.0 mA, in 
the sense of Gratton & Ortolani minus Minniti et al. Therefore, the larger equivalent widths 
with the warmer surface temperature are the probable cause of the higher [Fe/H] value they 
derived. Therefore, the detailed comparison of each elemental abundance between this study 
and Gratton & Ortolani may not be meaningful. 

Minniti et al. (1993) derived an even higher temperature of the star, T^s = 4400 K, which 
is 300 K warmer than our temperature, and, consequently, their iron abundance of the star 
from Fe I lines is about 0.45 dex higher than our value, but, interestingly, very close to the 
value we obtain using only the Fe II lines. Minniti et al. (1993) obtained logg — 1.0 for the 
star, using the ionization equilibrium of iron lines, which is about 0.3 dex larger than those of 
Gratton & Ortolani (1989) and this study. This is very worrisome, since Minniti et al. (1993) 
used the same temperature determination method, the same model atmospheres as this study 
and similar oscillator strength values. Their instrumental resolution was comparable to ours. 
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and they employed the program WIDTH rather than MOOG, but both programs yield the 
same results, in our experience, at least when the atomic data are identical. The equivalent 
widths are also in good agreement, as noted earlier. In particular, the agreement in the 
equivalent widths of weak lines (i.e., the lines on the linear part of the curve of growth) 
appears to be excellent. For lines with log(WA/A) < —5.2, the difference in equivalent 
widths is —1.7 ± 1.4 mA (7 lines), in the sense of Minniti et al. minus this study. 

Instead of comparing our results to those of Minniti et al. (1993) directly, we have 
rederived the surface temperature using the equivalent widths measured by Minniti et al. 
(1993). In their Table 2"^, Minniti et al. (1993) presented Fc I and Fc II line widths of the 
star 160. Using their Fe I line widths with log(WA/A) < —5.2, we ran MOOG for this star 
again. At this time, we use a Kurucz model atmosphere with Teflf = 4400 K, log^f = 1.0, and 
■f^turb = 2.0 km s~^ as our initial input model. Note that these stellar parameters are derived 
by Minniti et al. for their analyses of the star as shown in Table 9. We show the results from 
our MOOG run in Figure 5a. We obtain [Fe/H] = —2.18 ± 0.29 (36 lines) and this value is 
very similar to that of Minniti et al., [Fe/H] = —2.11. As can be seen in the Figure, however, 
the scatter in our analysis using the line list of Minniti et al. is much larger than that Minniti 
et al. estimated (c < 0.07 dex per line) and we were not able to reproduce their results using 
their input data. We find it necessary to adopt a different temperature. Since the Fe I line 
at A 6353. 84A deviates far from the mean value of the relation, we excluded this hne and ran 
MOOG again with T^s = 4275 K and obtain [Fc/H] = -2.30 ± 0.22 dex (Figure 5b). For 
comparison, we show the results from MOOG run with our stellar parameters, Teff — 4100 
K, log^ = 0.7, and Vturh = 1-8 km s"^ in Figure 5c. We obtain [Fe/H] = -2.50 ± 0.13. In 
the Figure, the agreement between results from our equivalent widths and those of Minniti 
et al. for lines with % < 3 eV is good. For lines with x > 3 cV, the scatter in results from 
Minniti et al. is rather large. It is likely that the Fe I lines employed by Minniti et al. may 
have been contaminated, and as a result they would have obtained a higher iron abundance 
than we have. Their accidental inclusion of Fe I lines suffering from line blending made their 
measurements to meet the ionization equilibrium at high surface gravity log g = 1.0. Using 
our line width measurements listed in Table 2, the criterion of the ionization equilibrium of 
iron lines is not satisfied until the surface gravity becomes as low as log g — 0.0 and this 
surface gravity value is not only 1 dex smaller than that of Minniti et al., but is implausible. 

We conclude that we do not understand why we cannot obtain the same results as did 
Minniti et al. (1993) using their data, and we conclude that the agreement between their 
[Fe/H] result and ours is probably fortuitous. 



^The titles of column 2 and 3 in Table 2 of Minniti et al. (1993) should be switched. 
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5. DISCUSSION 

In the previous section, we discussed the rationale by which we have derived elemental 
abundance estimates. From the ionized iron lines, we obtain [Fe/H] = —2.16 ± 0.02 {a = 
0.04) for M68 (based on internal errors only). Element-to-iron ratios are matched using 
comparable ionization states (except for oxygen, given its very high ionization potential). 

We explore the abundances of a variety of elements relative to iron. In the interests 
of economy and of interest, we focus on comparisons between M68 and the comparably 
metal-poor globular cluster M15, based on the results of Sneden et al. (1997, 2000), and 
then discuss these in turn compared to other ensembles of field stars and clusters. We select 
M15 because it is one of the few metal-poor globular clusters with comparable data for 
the neutron capture elements lanthanum, barium, and europium, and partly because of the 
hypothesized common origin (Yoon & Lee 2002). 

5.1. Mixing or Primordial Variations? O, Na, Mg, Al 

Many globular clusters appear to show anticorrelations between the abundances of oxy- 
gen and sodium, and of magnesium and aluminum. The subject was reviewed by Kraft 
(1994), and has been revisited by numerous authors. The approach often taken has been 
that these anticorrelations arise from deep mixing, whereby material whose chemical com- 
positions have been altered by proton captures within the CNO cycle have been brought to 
the stellar surface. This concept has become less plausible with the discovery that such anti- 
correlations are seen in relatively unevolved stars in the metal-poor clusters NGC 6397 and 
NGC 6752 (Gratton et al. 2001) and in the metal-rich cluster 47 Tuc (Carretta et al. 2004). 
The variations seen in some of these elements in Table 5 nonetheless warrant a quick look at 
these possible anti-correlations in M68, and we show them here in Figure 6, in comparison 
with the results for M15 from Sneden et al. (1997). 

There are three relatively interesting results here. First is that M68 does show a varia- 
tion, especially in sodium. The second is that the large range in sodium abundances is not 
matched by a correlated (or anti-correlated) variation in the oxygen abundances, as appears 
to be the case for M15. Finally, the post-AGB star 117, indicated by an open square in the 
Figure, is consistent with no signs of deep mixing, despite it clearly being the most evolved 
star in our study. The lack of an anti-correlation and the lack of variations seen in star 
117 are consistent with a primordial variation for the differences in elemental abundances, 
presumably due to pollution during the earliest stages of the cluster formation and evolution 
by AGB nucleosynthesis and mass loss. 
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5.2. The Other Light Elements: Si, Ca, and Ti Abundances 

Table 5 shows that the other hght elements, silicon, calcium, and titanium, do not show 
any detectable variation in their abundances. These elements are of considerable interest 
nonetheless because they provide us with an opportunity to compare the nucleosynthesis 
histories of clusters, and, in principle, a means to compare the rate of star formation and, 
possibly, relative ages. 

If star formation began everywhere in the Galaxy at about the same time, then the 
abundances of elements that emerge from Type II supernovae, including these hght elements 
and r-process elements such as europium, will be enhanced relative to the abundances of 
elements that emerge from nucleosynthesis sites that appear more gradually, including s- 
process elements from AGB stars, and iron-peak elements from Type la supernovae. The 
details are complex, but this basic picture is consistent with the basic behavior of these "a" 
elements' abundances (see the discussions of Wyse & Gilmore 1988; Wheeler et al. 1989, and 
Carney 1996). If we compute unweighted averages for silicon, calcium, and titanium, we find 
for M68 that [a/Fe] = +0.34 ±0.22, which is roughly consistent with that found in M15 and 
other globular clusters. In Table 10, we summarize the [Si/Fe], [Ca/Fe], [Ti/Fe], and the 
mean [a/Fe] ratios of M68, other globular clusters (Gratton 1987; Gratton & Ortolani 1989; 
Kraft et al. 1992, 1995, 1997, 1998; McWiUiam et al. 1992; Sneden et al. 1994, 1997, 2000b; 
Brown et al. 1997, 1999; Ivans et al. 1999; Shetrone & Keane 2000) and field stars with 
—2.50 < [Fe/H] < —1.90 studied by Fulbright (2000). (The errors are those of the mean.) 
In the Table, old inner halo clusters denote metal-poor inner halo clusters NGC 6287 ([Fe/H] 
= -2.01), NGC 6293 (-1.99), and NGC 6541 (-1.76) and their Galactocentric distances 
are 1.6 kpc, 1.4 kpc, and 2.2 kpc, respectively (Lee & Carney 2002). The mean values of 
[a/Fe] do not appear to vary much, and are probably consistent within the internal errors 
(the simple averages) and any lingering systematic errors. 

A closer examination of the abundances of these elements does reveal some interesting 
differences, however. The silicon abundance of M68 appears to be similar to that found in 
M15, and both are enhanced relative to those of other clusters and field stars. 

The calcium abundance of M68 is likewise similar to those found in M15, and both are 
similar to other halo stars and clusters. There may, however, be a modest difference 0.15 
dex) between M68, M15, and other "old halo" clusters and "younger halo" and thick disk 
clusters. Note that this contradicts the assignment of M68 to the "younger halo" category. 

On the other hand, the titanium abundance of M68 appears to be much lower than all 
other clusters. Is this effect real? The titanium abundance of the metal-poor RGB stars 
using the neutral titanium lines may suffer from NLTE effects, such as an over-ionization. 
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and the resultant Ti abundance may be spurious. However, our Ti abundance analyses using 
the Ti II lines also yield a lower titanium abundance for our program stars, indicating that 
they are truly titanium deficient. In Figure 7, we show comparisons of observed spectrum of 
the star 93 with those of synthetic spectra near Si I A 5665.56A and Ti I A 6258.11, 6258.71, 
and 626 1.11 A. In the Figure, the observed spectrum is presented by histograms, synthetic 
spectra with [Si/Fe] = +0.68, [Ti/Fe] = +0.06 by solid lines (see also Table 5), and synthetic 
spectra with [Si/Fe] = +0.30, [Ti/Fe] = +0.30, which elemental abundances have long been 
thought to be normal (if there exists normal elemental abundances among globular clusters) 
for globular clusters, by dotted fines. The observed spectrum clearly shows that the Si I 
absorption line is too strong to be [Si/Fe] = +0.30, while the Ti I absorption fines are too 
weak to be [Ti/Fe] = +0.30. So the deficiency appears to be real. 

While we have included this element because its abundances often track those of the 
lighter elements, including oxygen, magnesium, silicon, and calcium, titanium may also be 
considered to be an iron-leak element. Explosive nucleosynthesis calculations of the massive 
stars (Woosley & Weaver 1995) predict that one of the major sources of the SNe II titanium 
yield is ^^Cr via the consecutive electron capture processes. Further, their models predicted 
that SNe II with masses in the range 25 - 40 Mq are likely to overproduce Si compared with 
Ti. We turn, therefore, to a more detailed comparison between M68 and M15, concentrating 
on the light element silicon, thought to be mostly produced in Type II supernovae, and 
nickel, mostly associated with Type la supernovae. Sneden et al. (1997) also found a high 
silicon abundance in M15. Restricting the sample to only those stars with well-determined 
abundances, they found < [Si/Fe] > = +0.62 ± 0.06. Using again only the stars with well- 
determined titanium abundances, and comparing [Ti I/Fe I] and [Ti II/Fe II], the four stars 
studied by Sneden et al. (1997, 2000) resufied in <[Ti/Fe]> = +0.27 + 0.08 (the errors are 
all given here as errors of the mean). Thus for M15, [Si/Ti] is +0.35, whereas for M68 it is 
much higher, +0.76. For the iron-peak element nickel, our seven stars reveal <[Ni/Fc]> = 
—0.11 ± 0.03 [a = 0.08) dex. We are uncertain how to compare our results to those from 
Sneden et al. (1997, 2000), however. Sneden et al. (1997) found <[Ni/Fe]> = +0.14 + 0.03, 
based on 12 stars, which is very different than what we have found for M68. But Sneden 
et al. (2000) found <[Ni/Fe]> = -0.21 ± 0.02 {a = 0.04) dex, based on three of the same 
stars studied earlier. Sneden et al. (2000) commented on the difficulties in establishing truly 
reliable stellar parameters and abundances, and, alas, this difference is yet another aspect 
of that problem. 

We conclude our discussion of these elements with the possibifity that the stars we have 
studied in M68 sampled a relatively high end of the initial mass function and the resultant 
supernovae. 
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5.3. The Heavy Neutron Capture Elements Ba, La, and Eu 

The elements heavier than the iron-peak elements can not be efficiently produced by 
the charged-particle interactions due to the large Coulomb repulsion between the nuclei, 
and they are thought to be produced through both slow (s-) and rapid (r-) neutron cap- 
ture processes. The s-process occurs mainly in low- (1-3 M©) or intermediate-mass (4-7 
Mq) AGB stars, while the r-process is thought to occur in type II supernovae explosions. 
Therefore, comparisons between the r-process (europium) and s-process elements (barium 
and lanthanum) provide a clue to the history of the Galactic nucleosynthesis, since r- and s- 
processes are thought to occur in stars with very different masses and, therefore, in different 
evolutionary timescales. 

Burris et al. (2000) studied abundances of neutron capture elements in a large sample 
of metal-poor giants, finding that a large star-to-star variations in the neutron capture 
elemental abundances (see also Gilroy et al. 1988, McWilliam et al. 1995). They suggested 
that this scatter in neutron capture elemental abundances results from inhomogeneity of the 
proto-stellar material polluted by SNe II nucleosynthesis ejecta at early stage in the Galaxy's 
history. Whether this is the correct interpretation or whether the analyses themselves are 
partly responsible for the scatter remains to be seen (see Cayrel et al. 2004). In Figure 8, we 
show Ba, La, and Eu abundances for globular clusters (Brown et al. 1997, 1999; Gratton et 
al. 1986; Gratton 1987; Gratton & Ortolani 1989; Ivans et al. 1999, 2001; Kraft et al. 1998; 
Lee & Carney 2002; McWilliam et al. 1992; Shetrone & Keane 2000; Sneden et al. 1997, 
2000a, 2000b, 2004) and field stars (Burris et al. 2000) as a function of metallicity. Table 11 
summarizes mean Ba, La, Eu abundances for M68 and M15 (Sneden et al. 1997, 2000a), 
other globular clusters, and field stars with —2.50 < [Fe/H] < —1.90 (Burris et al. 2000). 

The comparison with M15 and other clusters and stars is not quite so simple, however, 
as it appears. Sneden et al. (1997) discovered that whereas the abundances of the lighter 
elements silicon, calcium, and titanium did not vary from star to star, those of the neutron 
capture elements did, and by factors of four to five. Like the variation in [Si/Ti] seen 
between M68 and other clusters, the variations in [Ba/Fe] and [Eu/Fe] seen by Sneden et al. 
(1997, 2000) warn us that we cannot assume all supernovae produce the same heavy element 
abundance yields. In the case of M15, the lack of variations in the abundances of the fight 
elements would predict in such a simple model that barium and other s-process elements 
and the r-process element europium should likewise be constant, contrary to what was seen. 

How does M68 fit into this picture? In Figure 9 we reproduce the [Ba/Fe] vs. [Eu/Fe] 
abundances found by Sneden et al. (1997). The dotted line shows the solar abundances, and 
the dashed fine, displaced from the solar relation by 0.41 dex, shows the mean behavior of 
the stars in M15, which are plotted as dots. As Sneden et al. (1997) noted, the behavior of 
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the stars in M15 is consistent with an enhanced r-process contribution to the abundances 
of both barium and europium in M15. That process produced different absolute amounts of 
the neutron capture elements in stars in Ml 5, but the process itself yielded the same relative 
abundances of the two elements. Our results for the stars in M68 are consistent with this 
general picture, in that the barium and europium abundance ratio is constant, and essentially 
the same in M68 and in M15. Sneden et al. (1997) found that the "low barium" and "high 
barium" groups of stars in M15 had similar [Ba/Eu] ratios, —0.41 ± 0.03 (cr = 0.13) dex for 
all eighteen stars. Our seven stars, including the post-AGB star 117, yield < [Ba/Eu] > = 
— 0.43±0.05 ((7 = 0.13) dex. What is different between the two clusters, however, is that both 
barium and europium are lower in abundance in M68 than the stars in M15, even the "low 
barium" stars. Thus M68 was slightly less enriched in r-process nucleosynthesis, relative to 
iron. We could interpret our results in an alternative fashion, that iron is enhanced relative 
to the r-process elements, in this case including barium and europium. 

For [La/Eu] we must turn to the smaller sample of Sneden et al. (2000), who studied 
three stars in M15 with very different [Ba/Fe] ratios, —0.24, +0.05, and +0.05. The [La/Fe] 
values also vary greatly for these three stars: —0.02, +0.38, and +0.61, yet the variation 
seen in [La/Eu] is much smaller, and consistent with the measurement uncertainties: —0.33, 
—0.45, and —0.26, for a mean value of —0.35 ± 0.06 (a = 0.10) dex. We were able to 
estimate [La/Eu] values for three of our program stars, and we find < [La/Eu] > = —0.38 ± 
0.10 ((7 = 0.17) dex. Thus for both lanthanum and barium, it appears the r process is 
dominating the nucleosynthesis enrichment in both M68 and in M15, and yielding constant 
relative abundances of the elements. But as in the case of barium and europium, there is an 
underabundance of lanthanum relative to iron in M68 compared to M15. 



5.4. Summary of Differences Between M68 and M15 and Other Clusters 

In Figure 10, we summarize graphically the many similarities between the abundances of 
various elements, relative to iron, in M68 compared to M15. Perhaps the two most striking 
differences are in the much lower titanium abundances in M68, which, as we have seen, 
might be explained by a greater contribution from more massive stars and their supernovae 
events in M68 relative to M15. On the other hand, while it appears that the r-process has 
dominated the production of the "traditional" s-process elements lanthanum and barium, 
and that the [La/Eu] and [Ba/Eu] ratios are the same for all the stars in both clusters, 
the overall r-process enrichment varies within M15 (Sneden et al. 1997) and between M68 
and M15. The clusters have clearly experienced somewhat different chemical enrichment 
histories, and it will be interesting to see if models of supernovae enrichment ultimately 
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prove successful in explaining these differences and what this means for the histories and, 
perhaps, ages of these clusters. 

Finally, in Figure 11 we show the results for individual stars in the abundance ratio of 
barium and europium, plotted against [Ca/H] rather than [Fe/H]. As Sneden et al. (1997) 
commented, calcium is thought to be produed in Type II supernovae, compared to iron, 
whose abundance is more strongly affected by contributions from Type la supernovae. Also 
following Sneden et al. (1997), we show the levels expected from "pure" r-process nucle- 
osynthesis, "pure" s-process nucleosynthesis, and the intermediate case found to exist in the 
solar system. At lower [Ca/H] levels, [Ba/Eu] ratios more consistent with r-process domina- 
tion are apparent, which is not terribly surprising. But M68 appears to be one of the most 
extreme cases. 



6. CONCLUSIONS 

A chemical abundance study of seven giant stars in M68 has been presented. We 
estimate the stars' temperatures using the abundances derived from Fe 1 lines with differing 
excitation potentials, but the gravity has been derived using available photometry rather 
than a comparison between abundances from pressure-insensitive Fe I lines and pressure- 
sensitive Fe II hue. The "spectroscopic" gravities do not agree with the photometric gravities, 
suggesting non-LTE effects are present. Using only the Fe II hues, which should be less 
sensitive to non-LTE, and photometric gravities, we find [Fe/H] = —2.16 ± 0.02 {a = 0.04. 
We have compared our results to those of Gratton & Ortolani (1989), who found a higher 
value. We attribute the difference to the lower-resolution spectra available to them at the 
time. Our results do agree well with those obtained by Minniti et al. (1993), but we regard 
the agreement as accidental since we could not reproduce their results using their data. 

We determine element-to-iron ratios using neutral vs. neutral and ionized vs. ionized 
lines to again minimize non-LTE effects. We find a large range in sodium abundances, 
but no significant range in oxygen abundances. Further, the post-AGB star M68-117 does 
not appear to show any enhancement of sodium. These results are not consistent with 
deep mixing being the cause of the variations among the light elements oxygen, sodium, 
magnesium, and aluminum. 

There are two notable differences between M68 and the comparably metal-poor cluster 
Ml 5. While both show enhanced [Si/Fe] ratios relative to other clusters and comparably 
metal-poor field stars, M68 is quite deficient in titanium compared to M15 or any other 
cluster. It is possible that this arose because nucleosynthesis enrichment of the stars in M68 
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was provided by supernovae resulting from the deaths of somewhat more massive progenitors. 
This would be difficult to reconcile with the nominal younger age for M68 compared to M15, 
based on the morphologies of their horizontal branch stars' distribution in temperature/color. 
Perhaps age is not the only "second parameter"? 

The second interesting difference is that [La/Eu] and [Ba/Eu] ratios are similar for the 
stars in M68 and in M15 (and in both its "high barium" and "low barium" groups of stars — 
Sneden et al. 1997). This suggests that in both clusters the r process is making a major 
contribution to the abundance levels of the "traditional" s-process elements lanthanum and 
barium. However, whatever the process is, it is not contributing quite as much to the neutron 
capture elements' abundances in M68 as in M15. 
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Table 1. Journal of observations. 
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i"W" from Walker (1994). "A" from Alcaino (1977). T' from Harris (1975). 
^Positional data from 2MASS (Cutri et al. 2000). 
^Walker (1994). 

^Frogel, Persson, & Cohen (1983). 

52MASS photometric data using the CIT system (Cutri et al. 2000). 
^An ultraviolet bright star (Zinn, Newell, & Gibson 1972). 



GIANT STARS IN M68 



24 



Table 2. Equivalent widths. 
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Table 2 — Continued 
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Table 2 — Continued 
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Table 2 — Continued 
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62 


63 


5 


6141.73 


Ba II 


0.700 


-0.080 


131 


109 


120 


113 


117 


113 


103 


5 


6496.91 


Ba II 


0.600 


-0.380 


127 


109 


119 


122 


119 


113 


98 


5 


6390.48 


La II 


0.321 


-1.410 


5 




5 




6 






4 


6774.27 


La II 


0.126 


-1.750 


7 








6 






4 


6645.06 


Eu II 


1.380 


0.204 


6 


5 


8 


9 


7 


10 




22 


7217.55 


Eu II 


1.230 


-0.301 


3 




7 


8 








22 


7426.57 


Eu II 


1.279 


-0.149 


5 














22 
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Table 3. References for log gf values. 



Kef. No. 


Reference 


tilem. 


1 


Lambert (1978) 


\0 11 

J 


2 


Ivans et al f2001) 


Na I 


3 


Kraft ct al. (1992) 


Na I 


4 


Sneden et al. (2004) 


Mg 1, Si I, Ca I, 






Ti 1, Ni 1, La 11 


5 


Sneden et al. (1997) 


Al I, Ba II 


6 


Garz (197i) 


Si I 


7 


Smith fl981) 


Ca I 


8 


Blackwell et al. (1982a) 


Ti I 


9 


Blackwell et al. (1983) 


Ti I 


10 


Blackwell et al. (1986b) 


Ti I 


11 


Fuhr et al. (1988a) 


Ti II 


12 


Prochaska & McWilliam (2000) 


Sc I, Mn I 


13 


Blackwell et al. (1982b) 


Fe I 


14 


Blackwell et al. (1982c) 


Fe I 


15 


Blackwell et al. (1986a) 


Fe I 


16 


O'Brian et al. (1991) 


Fe I 


17 


Bard et al. (1991) 


Fe I 


18 


Bard & Kock (1994) 


Fe I 


19 


Biemont (1991) 


Fe II 


20 


Fuhr et al. (1988b) 


Ni I 


21 


Kurucz (1993) 


Cu I 


22 


Biemont et al. (1982) 


Eu II 
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Table 4. Model atmosphere parameters. 



Star 






rcir(K) 






^ogg 
(phot) 


^ogg 
(spec) 


I'turb 

(km s~^) 


{B - vy 


{V 

(FPC83) 


-Ky 

(2MASS) 


(FPC83)2 


(spec) 


93 


4220 


4215 


4240 


4287 


4200 


0.7 


0.1 


1.95 


117 


4370 


4340 


4335 


4501 


4300 


0.8 


0.3 


1.80 


160 


4200 


4200 


4200 


4329 


4100 


0.7 


0.0 


1.80 


89 


4200 


4200 


4215 


4279 


4175 


0.7 


0.0 


1.90 


412 


4275 


4290 


4300 


4372 


4200 


0.8 


0.1 


1.75 


350 


4390 




4340 




4300 


1.0 


0.3 


1.60 


48 


4450 




4450 




4325 


1.1 


0.3 


1.75 



^Using relations given by Alonso et al. (1999). 
^Progelet al. (1983). 
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Table 5. Elemental abundances of M68 with photometric logg'. 





93 


117 


160 


89 


412 


350 


18 


Mean 




[te/ n\i 


— z.ou 


— Z.Oo 


— Z.OD 


9 c^l 
— Z.Oi 


9 c^S 
— Z.Oo 


9 K'i 
— Z.Oo 


9 fi7 
— Z.D / 


— Z.Oo 


— Z.OD 


n 




oo 


49 
^Z 


^O 


49 
^Z 


OO 


OO 


7 


fi 
D 


(7 


u.uo 


U.UD 


u.uo 


n nc; 
u.uo 


n OR 

U.UD 


u.uo 


OR 
U.UD 


OQ 

u.uy 


Ofi 
U.UD 


[te/ tijn 


— Z.l / 


9 1 R 


9 19 

— z.iz 


9 14 
— Z. 141 


9 1 fi 
— Z.ID 


9 1 
— Z.IO 


9 9"? 
— z.zo 


9 1 fi 
— Z.ID 


9 1 fi 
— Z.ID 


n 
11 


Q 
o 


r. 


o 




Q 
o 


Q 

o 


Q 
o 


7 


fi 


(7 


U.UO 


u.uo 


n 09 
u.uz 


n nfi 
u.uo 


n n^; 
u.uo 


09 
U.UZ 


Ofi 
U.UD 


0"? 

u.uo 


04 
U.U4 


rn /FpI 
[u/rej 






-I-n 9'^ 


4-0 4'? 


4-0 '^7 


4-0 ^1 


4-0 'iQ 


4-0 49 


4-0 49 


n 


1 
1 


9 

z 


9 

z 


9 
z 


1 
1 


1 
1 


9 
z 


7 


D 


(7 




n on 
u.zu 


U.Ui 


n 9n 
u.zu 






n nc; 
u.uo 


n 1 9 
u. iz 


n 1 q 
u. lo 


[No /"Pel 


4-0 fi9 


n n7 

— u.u / 


n 9"^ 

— U.ZO 


n 97 

— U.Z / 


4-0 "^S 


4-0 

T^U.OO 


4-0 "^Q 
T^u.oy 


4-0 01 
T^U.ZO 


4-0 94 


n 


Q 
O 


1 
1 


9 

z 


1 
1 


Q 
O 


9 
Z 


1 
i 


7 


R 
D 


(7 


u.u / 




u.uu 




n na 

U.Uo 


O Oc; 

u.uo 




O /lO 
U.4U 




u.oy 


fMo- /Pel 

[ivig/rej 


4-0 1 9 


-i-fl "^Q 
T^u.oy 


T^U.OU 


4-0 47 


4-0 'i'^ 


4-0 1Q 


4-0 4'i 


4-0 'i'^ 
T^U.OO 


4-0 'i'^ 
T^U.OO 


n 


z 


z 


9 

z 


9 

Z 


9 
Z 


1 


9 
Z 


7 
( 


fi 
D 




u.uo 


U.Ul 


n ni 

U.Ul 


n ni 

U.Ul 


n 1 n 

U. iU 




n nQ 
u.uo 


n 1 /I 

U. 14 


n 1 fi 

U. ID 


r A 1 /PpI 


-1-1 1 s 








4-1 nn 

T^l.UU 


4-1 07 

T^l.U / 




4-1 ns 

T^l.Uo 


4-1 ns 


n 


z 








9 

z 


9 
Z 




q 
o 


q 
o 


(7 


u.uo 








07 
U.U / 


07 
U.U ( 




OQ 

U.uy 


OQ 

U.uy 


fQi /"Pol 


4-0 fiS 
T^U.Oo 


-i-fi fill 


4-0 77 




4-0 7^ 


4-0 79 




4-0 71 


4-0 7*^ 


n 


q 
o 




z 


1 
i 




Q 
O 


1 
1 







A 


(7 


0.05 


0.04 






0.05 






0.05 


0.04 


[Ca/Fe] 


+0.33 


+0.27 


+0.29 


+0.37 


+0.30 


+0.23 


+0.32 


+0.30 


+0.31 


n 


9 


10 


10 


7 


7 


9 


5 


7 


6 


a 


0.05 


0.07 


0.04 


0.08 


0.04 


0.07 


0.05 


0.04 


0.05 


[Sc/Fe]n 


+0.14 


+0.19 


+0.16 


+0.14 


+0.13 


+0.07 


+0.13 


+0.14 


+0.13 


n 


4 


4 


3 


4 


4 


4 


4 


7 


6 


(7 


0.14 


0.20 


0.25 


0.20 


0.18 


0.18 


0.18 


0.04 


0.03 


[Ti/Fe]i 


+0.06 


+0.07 


+0.09 


+0.03 


-0.02 


-0.02 


+0.01 


+0.03 


+0.02 


n 


8 


6 


7 


8 


6 


5 


4 


7 


6 


(7 


0.07 


0.05 


0.07 


0.05 


0.04 


0.06 


0.05 


0.04 


0.05 
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Table 5 — Continued 





93 


117 


160 


89 


112 


350 


18 


Mean 




[Ti/Fe]„ 


-0.12 


+0.09 


+0.01 


-0.15 


-0.06 


-0.08 




-0.05 


-0.08 


n 


2 


2 


2 


2 


2 


2 




6 


5 


(7 


0.03 


0.12 


0.08 


0.05 


0.01 


0.02 




0.09 


0.06 


[Ti/Fe]„,ean 


-0.03 


+0.08 


+0.05 


-0.06 


-0.04 


-0.05 


+0.01 


-0.01 


-0.03 


[Mn/Fe] 


-0.25 


-0.22 


-0.30 


-0.27 


-0.29 


-0.26 


-0.33 


-0.27 


-0.28 


n 


1 


1 


1 


1 


1 


1 


1 


7 


6 


cr 
















0.04 


0.03 


[Ni/Fe] 


-0.04 


-0.21 


-0.01 


-0.09 


-0.09 


-0.22 


-0.14 


-0.11 


-0.10 


n 


4 


3 


3 


3 


3 


3 


2 


7 


6 


a 


0.11 


0.11 


0.08 


0.09 


0.09 


0.09 


0.00 


0.08 


0.07 


[Cu/Fe] 


-0.74 




-0.83 


-0.86 




-0.85 




-0.82 


-0.82 


n 


1 




1 


1 




1 




4 


4 


(7 
















0.05 


0.05 


[Ba/Fe]„ 


-0.28 


-0.46 


-0.39 


-0.40 


-0.30 


-0.19 


-0.32 


-0.33 


-0.31 


n 


3 


3 


3 


3 


3 


3 


3 


7 


6 


a 


0.09 


0.12 


0.06 


0.10 


0.12 


0.14 


0.04 


0.09 


0.08 


[La/Fe]„ 


-0.23 




-0.40 




-0.19 






-0.27 


-0.27 


n 


2 




1 




2 






3 


3 


a 


0.14 








0.03 






0.11 


0.11 


[Eu/Fe]„ 


+0.06 


+0.01 


+0.18 


+0.28 


+0.11 


+0.37 




+0.17 


+0.20 


n 


3 


1 


2 


2 


1 


1 




6 


5 


(7 


0.08 




0.13 


0.12 








0.14 


0.13 



^Without the star 117. 
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Table 6. Elemental abundances of M68 spectroscopic log g. 



93 



117 



160 



89 



112 



350 



18 



Mean Meaiii 



[Fe/H]i 
n 

(7 

[Fe/H]„ 
n 
a 



-2.40 
43 
0.05 
-2.37 
3 

0.02 
-2.39 



-2.30 
38 
0.06 
-2.33 
5 

0.05 
-2.32 



-2.46 
42 
0.05 
-2.35 
3 

0.02 
-2.41 



-2.40 
45 
0.06 
-2.36 

4 

0.06 
-2.38 



-2.47 
42 
0.06 
-2.39 
3 

0.04 
-2.43 



-2.43 
35 
0.05 
-2.39 
3 

0.02 
-2.41 



-2.56 
35 
0.07 
-2.52 
3 

0.06 
-2.54 



-2.43 
7 

0.08 
-2.39 
7 

0.06 
-2.43 



-2.45 
6 

0.06 
-2.40 
6 

0.06 
-2.41 



[0/Fe] 
n 

(7 

[Na/Fe] 
n 

cr 

[Mg/Fe] 
n 
cr 

[Al/Fe] 
n 

(7 

[Si/Fe] 
n 

[Ca/Fe] 
n 
cr 

[Sc/Fe]„ 
n 
cr 

[Ti/Fe]i 
n 



+0.38 
1 

+0.65 
3 

0.08 
+0.19 
2 

0.03 
+1.16 
2 

0.06 
+0.62 
3 

0.06 
+0.36 
9 

0.04 
+0.23 
4 

0.19 
+0.06 



+0.43 
2 

0.20 
-0.04 
1 

+0.48 

2 

0.06 



+0.63 
2 

0.04 
+0.31 
10 
0.08 
+0.22 
4 

0.22 
+0.08 
6 



+0.23 
2 

0.01 
-0.21 
2 

0.00 
+0.60 

2 

0.08 



+0.67 
1 

+0.31 
10 
0.07 
+0.33 
3 

0.34 
+0.06 
7 



+0.42 
2 

0.21 
-0.24 
1 

+0.60 

2 

0.08 



+0.39 
7 

0.09 
+0.23 
4 

0.26 
+0.08 



+0.39 
1 

+0.38 
3 

0.07 
+0.43 

2 

0.21 
+0.98 
2 

0.07 
+0.66 
3 

0.06 
+0.31 
7 

0.06 
+0.24 
4 

0.24 
-0.04 
6 



+0.51 
1 

+0.56 
2 

0.05 
+0.32 
1 

+1.04 
2 

0.07 
+0.64 
1 

+0.24 
9 

0.07 
+0.14 
4 

0.22 
-0.02 
5 



+0.60 
2 

0.06 
+0.40 
1 

+0.53 
2 

0.13 



+0.32 

5 

0.05 
+0.20 
4 

0.22 
+0.01 

4 



+0.42 
7 

0.11 
+0.21 
7 

0.37 
+0.45 
7 

0.15 
+1.06 
3 

0.09 
+0.64 
5 

0.02 
+0.32 
7 

0.05 
+0.23 
7 

0.06 
+0.03 
7 



+0.42 
6 

0.13 
+0.26 
6 

0.39 
+0.45 

6 

0.17 
+1.06 
3 

0.09 
+0.65 
4 

0.02 
+0.32 
6 

0.05 
+0.23 
6 

0.06 
+0.03 
6 
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Table 6 — Continued 





93 


117 


160 


89 


412 


350 


48 


Mean 


Mean^ 




0.07 


0.05 


0.09 


0.06 


0.04 


0.06 


0.06 


0.05 


0.05 


[Ti/FeliT 


—0.10 


+0.08 


+0.06 


—0.15 


—0.02 


—0.06 




—0.03 


—0.05 


n 


2 


2 


2 


2 


2 


2 




6 




(J 


n 02 


13 


10 


06 


00 


03 




09 


08 

W . WtJ 


[Ti/Fel 


—0 02 


-1-0 08 


-1-0 06 


—0 04 


— 03 


—0 04 

w . w^ 


+0 01 


+0 00 


—0 01 

W . W X 


fMn/Fel 


-0.24 


-0.19 


-0.31 


-0.25 


-0.30 


-0.26 


-0.33 


-0.27 


-0.28 




1 


1 


1 


1 


1 


1 


1 


7 

f 


w 


















05 


04 

w . W^ 


[Ni/Fel 


-0.04 


-0.22 


-0.13 


-0.16 


-0.18 


-0.27 


-0.19 


-0.17 


-0.16 




4 


3 


3 


3 


3 


3 


2 


7 

f 


R 
w 


0" 


0.11 


0.11 


0.09 


0.09 


0.09 


0.09 


0.00 


0.07 


0.08 


[Cu/Fe] 


-0.76 




-0.90 


-0.88 




-0.87 




-0.86 


-0.86 


n 


1 




1 


1 




1 




4 


4 


cr 
















0.06 


0.06 


[Ba/Fe]n 


-0.18 


-0.44 


-0.24 


-0.31 


-0.18 


-0.12 


-0.24 


-0.24 


-0.21 


n 


3 


3 


3 


3 


3 


3 


3 


7 


6 


(7 


0.11 


0.11 


0.05 


0.10 


0.12 


0.13 


0.06 


0.11 


0.07 


[La/Fe]„ 


-0.20 




-0.33 




-0.15 






-0.23 


-0.23 


n 


2 




1 




2 






3 


3 


cr 


0.13 








0.02 






0.10 


0.10 


[Eu/Fe]„ 


+0.06 


+0.01 


+0.24 


+0.28 


+0.15 


+0.39 




+0.18 


+0.22 


n 


3 


1 


2 


2 


1 


1 




6 


5 


a 


0.08 




0.12 


0.11 








0.14 


0.13 



1 Without the star 117. 
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Table 7. Differences in elemental abundances using photometric and spectroscopic surface 

gravities. 



Elem. 


[el/ Fe] phot - [el/ Fe] spec 


rrp^ /Til 

[te/ti] 


+U.ZO ± U.Uo 


[0/Fe] 


+0.01 ± 0.01 


[Na/Fe] 


-0.02 ± 0.01 


[Mg/Fe] 


-0.10 ± 0.02 


[Al/Fe] 


+0.02 ± 0.01 


[Si/Fe] 


+0.07 ± 0.03 


[Ca/Fe] 


-0.02 ± 0.01 


[Sc/Fe]„ 


-0.09 ± 0.04 


[Ti/Fe]i 


+0.00 ± 0.03 


[Ti/Fe]ii 


-0.02 ± 0.02 


[Ti/Fe]iYiean 


-0.01 ± 0.01 


[Mn/Fe] 


-0.01 ± 0.02 


[Ni/Fe] 


+0.05 ± 0.04 


[Cu/Fe] 


+0.03 ± 0.03 


lBa/Fe]ii 


-0.10 ± 0.04 


[La/Fe]„ 


-0.05 ± 0.02 


[Eu/Fe]„ 


-0.02 ± 0.03 
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Table 8. Abundance dependencies on model atmosphere. 



Elem. STes Slogg Svturh 



±80 ± 0.3 ± 0.2 

(K) (km s-^) 



[Fe/H]i 


± 


0.12 


T 


0.03 


T 


0.03 


[Fe/H]„ 




0.06 


± 


0.10 


± 


0.02 


[0/Fe] 


± 


0.04 


± 


0.05 




0.03 


[Na/Fe] 


± 


0.03 


T 


0.03 


± 


0.03 


[Mg/Fe] 


T 


0.05 


T 


0.05 


± 


0.01 


[Al/Fe] 


T 


0.06 


T 


0.01 


± 


0.03 


[Si/Fe] 


T 


0.10 


± 


0.01 


± 


0.04 


[Ca/Fe] 




0.04 




0.04 


± 


0.01 


[Sc/Fe]n 


± 


0.06 


T 


0.02 




0.01 


[Ti/Fe]i 


± 


0.03 


T 


0.02 


± 


0.03 


[Ti/Fe]„ 


± 


0.04 


± 


0.00 


± 


0.02 


[Mn/Fe] 


T 


0.02 




0.03 


± 


0.04 


[Ni/Fe] 




0.02 


± 


0.01 


± 


0.01 


[Cu/Fe] 


± 


0.01 


± 


0.01 


± 


0.02 


[Ba/Fe]„ 


± 


0.08 




0.02 




0.09 


[La/Fe]„ 


± 


0.09 




0.01 


± 


0.02 


[Eu/Fe]„ 


± 


0.06 


± 


0.00 


± 


0.02 
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Table 9. Stellar parameters and elemental abundances of the star 160. 



This study Gratton & Ortolani Minniti et al. 

(1989) (1993 & 1996) 





4100 


4329 


4400 


(K) 










0.7 


0.75 


1.0 


method 


phot. 


phot. 


spec. 




1.8 


1.6 


2.0 


(km sec~^) 








[Fe/H]i 


-2.56 


-1.94 


-2.11 


[Fe/H]„ 


-2.17 




-2.14 


[0/Fe] 


+0.23 


+0.22 


+0.25 


[Na/Fe] 


-0.23 


-0.04 


-0.08 


[Mg/Fe] 


+0.50 


+0.07 




[Si/Fe] 


+0.77 


+0.41 




[Ca/Fe] 


+0.29 


+0.25 




[Ti/Fe] 


+0.05 


+0.43 




[Ni/Fe] 


-0.01 


-0.15 




[Ba/Fe] 


-0.39 


-0.33 
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Table 10. Comparisons of a-element abundances. 



[Si/Fc] [Ca/Fc] [Ti/Fc] [a/Fc] 



M68 


0, 


.73 


± 


0, 


.02 


0. 


,31 


± 


0. 


,02 


-0.03 ± 0.04 


0. 


,34 


± 


0, 


.22 


5,7,7 


MI52 


0, 


.62 


± 


0, 


.06 


0. 


,24 


± 


0. 


,01 


0.27 ± 0.08 


0. 


,38 


± 


0, 


.12 


11,18,4 


old halo 


0, 


.38 


± 


0, 


.05 


0. 


,30 


± 


0. 


,03 


0.33 ± 0.03 


0. 


,34 


± 


0, 


.02 


15 


old inner halo^ 


0, 


.56 


± 


0, 


.05 


0. 


,26 


± 


0. 


,07 


0.15 ± 0.06 


0. 


,32 


± 


0, 


.03 


3 


younger halo 


0, 


.29 


± 


0, 


.09 


0. 


,14 


± 


0. 


,06 


0.17 ± 0.07 


0. 


,20 


± 


0, 


.05 


7 


younger halo^ 


0, 


.38 


± 


0, 


.09 


0. 


,21 


± 


0. 


,04 


0.25 ± 0.06 


0. 


,28 


± 


0, 


.02 


5 


disk 


0, 


.39 


± 


0, 


.08 


0. 


,16 


± 


0. 


,07 


0.31 ± 0.08 


0. 


,29 


± 


0, 


.04 


4 


Field^ 


0, 


.52 


± 


0, 


.04 


0. 


,35 


± 


0. 


,03 


0.30 ± 0.04 


0. 


,39 


± 


0, 


.02 


9, 20, 20 



^Number of stars (M68; M15; Field) or clusters. 

^We used [Si/Fe] and [Ca/Fe] from Sneden et al. (1997). For [Ti/Fe], we also employed 
data from Sneden et al. (2000). We did not include any stars with uncertain measures 
(marked with a colon in their results. 

^NGC 6287, NGC 6293, & NGC 6541 (Lee & Carney 2002). 

^Without Palomar 12 and Rupcrcht 106. 

^Mean abundances for field stars with —2.50 < [Fe/H] < —1.90. 
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Table 11. Comparisons of neutron capture element abundances. 



[Ba/Fe] n [La/Fe] n [Eu/Fe] 



M68 


-0, 


.31 


± 


0, 


.08 




-0, 


.27 


± 


0.11 




+0, 


.20 


± 


0, 


.13 


7,3,6 


M15-alP 


+0, 


.10 


± 


0. 


.21 














+0 


.49 


± 


0, 


.20 


18,0,18 


MlS-high Ba=^ 


+0, 


.23 


± 


0. 


.09 














+0 


.66 


± 


0, 


.08 


10,0,10 


MI5-I0W Ba^ 


-0, 


.11 


± 


0, 


.06 














+0 


.29 


± 


0, 


.09 


8,0,8 


old halo 


+0, 


.11 


± 


0, 


,23 


16 


+0, 


.22 


± 


0.18 


5 


+0 


.45 


± 


0, 


.13 


8 


old inner halo^ 


+0, 


.21 


± 


0. 


,20 


3 


+0, 


.18 


± 


0.08 


3 


+0 


.39 


± 


0, 


.07 


3 


younger halo 


-0, 


.03 


± 


0, 


.33 


5 


+0, 


.17 


± 


0.12 


3 


+0 


.48 


± 


0, 


.35 


4 


younger halo^ 


+0, 


.05 


± 


0. 


,34 


3 


+0. 


09 






1 


+0, 


.56 


± 


0, 


.02 


2 


disk 


-0, 


.25 


± 


0, 


.22 


4 












+0. 


33 








1 


Field^ 


+0, 


.18 


± 


0. 


.11 


27 


+0, 


.11 


± 


0.08 


14 


+0, 


.39 


± 


0, 


.09 


17 



^Number of stars (M68; M15; Field) or clusters. 

^Sneden et al. (1997). 

3 "High Ba" group (Sneden et al. 1997). See Figure 9 
^"Low Ba" group (Sneden et al. 1997). See Figure 9 

^NGC 6287, NGC 6293, & NGC 6541 (Lee & Carney 2002). 
^Without Palomar 12 and Rupercht 106. 

^Mean abundances for field stars with -2.50 < [Fe/H] < -1.90 (Burris et al. 2000). 
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Fig. 1. — HB type versus [Fe/H] (HB isochrones) of "halo" globular clusters (Da Costa & 
Armandroff 1995) and M68. Crosses are "old halo" clusters and open circles are "younger 
halo" clusters, where large crosses and open circles denote clusters studied employing high- 
resolution spectroscopy (see also Figure 2). HB isochrones for At — 0.0, —1.1, and —2.2 Gyr 
(with respect to the mean age of "old halo" globular clusters) are also shown with dashed 
hues (Rey et al. 2001). 
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Fig. 2. — Color- magnitude diagram of M68 with photometry from Walker (1994), showing 
the positions of our program stars. 
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Fig. 3. — A comparison of equivalent widths of the star 160 with those of Gratton & Ortolani 
(1989) and Minniti et al (1993). 
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Fig. 4. — Elemental abundances studied in this work plotted against Teflf- Note that the 
ordinate of the top panel is [Fe/H], while the ordinates of the remaining panels are [el/Fe] 
ratios. Neutral atoms are presented by filled circles and singly ionized atoms by open circles. 
Dashed lines are the mean abundances of each element. The mean Fe II and Ti II abundances 
are presented by dotted lines. 
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Fig. 5. — (a) The log n(Fe) versus excitation potential relation of the star M68-160 (1260) 
using Fe I equivalent widths of Minniti ct al. (1993) with log(WA/A) < -5.2 for T^s = 4400 
K, \ogg = 1.0, and t'tmb = 2.0 km s^^. Crosses are for Fe I lines with Wa > 10 niA and open 
squares are for Wa < 10 mA from Minniti et al. (1993). The dotted and the dashed lines 
represent the least square fits to the data with and without the Fe I A 6538.84 A, respectively. 
Filled circles show the relation using Fe I equivalent widths of this study and the sohd hne 
is for the least square fit to the data, (b) Same as (a), but using the model atmosphere with 
Teff = 4275 K, log^ = 1.0, and ^;turb = 2.0 km s"^ and the Fe I A 6538.84 A hne is not used, 
(c) Same as (b), but using the model atmosphere with Teff = 4100 K, logg^ = 0.7, and t'turb 
= 1.8 km s~^. 
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Fig. 6. — A comparison of [Na/Fe] vs. [0/Fe] for M15 (dots) vs. those of M68 (circles). The 
post-AGB star M68-117 is marked with an open square. 



GIANT STARS IN M68 



45 



>^ 1.00 
•I— I 



0.95 - 



1 1 




s 

1 1 


i I 

, 1 , 



5664 



5666 

A(A) 




6258 



6260 



x(A) 



6262 



Fig. 7. — Comparisons of observed spectrum of the star M68-93 with those of synthetic 
spectra near Si I A 5665.56A and Ti I A 6258.11, 6258.71, and 6261. llA. The observed 
spectrum is presented by histograms, synthetic spectra with [Si/Fe] = +0.68, [Ti/Fe] = 
+0.06 by sohd hues (see also Table 5), and synthetic spectra with [Si/Fe] — +0.30, [Ti/Fe] 
= +0.30 by dotted lines. 
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Fig. 8. — Abundances of the neutron capture elements. Crosses are old halo clusters, open 
circles are younger halo clusters, open squares arc thick disk clusters, and grey dots are field 
stars (Burris et al. 2000). M68 is represented by filled circles. 
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Fig. 9. — The abundances of barium and europium relative to iron in M15 (dots) and 
M68 (circles). The post-AGB star M68-117 is marked by an open square. The dotted 
line represents the solar abundance pattern. The dashed line is the approximate mean for 
the stars in M15. Despite varying levels of neutron capture abundances in the cluster, the 
[Ba/Eu] ratio appears to be constant. That is true as well for M68. 
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Fig. 10.— 



Comparison of elemental abundances between M68 and M15. 
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Fig. 11. — Comparisons of [Ba/Eu] ratios as a function of [Ca/H] of our program stars in 
M68 and individual stars in other globular clusters. 



